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I. INTRODUCTION 


The thermal radiation of photons and dileptons from a quark-glnon plasma (QGP) is 
one of the most promising signatures for the formation of a QGP in relativistic heavy ion 
collisions [Q]. Whereas real photons can be produced to lowest order aag only via the 
participation of a gluon, dileptons are produced to lowest order o? by quark-antiquark 
annihilation into a virtual photon (Born term) However, in the case of a small invariant 
photon mass where E is the energy and p = |p| the momentum of the photon, 

radiative corrections of the order a^ag become increasingly important. These corrections 
have been considered using perturbative QGD at finite temperature [@,|1. In contrast to the 
production of real photons the dilepton production rate turns out to be infrared finite in the 
case of a vanishing quark mass. There is a cancellation of the infrared singularities of real 
and virtual contributions, the latter appearing only for the production of virtual photons 

However, also medium effects leading to an effective quark mass m* by the interaction of 
the quark with the heat bath can screen these infrared divergences. Altherr and Ruuskanen 
0 suggested that the effective quark mass simply replaces the invariant photon mass as 
infrared cutoff if M < m*. 

Medium effects can be included consistently using the resummation technique of Braaten 
and Pisarski [^. This method has been applied to a number of interesting quantities of the 
QGP leading to gauge invariant and infrared finite results that are complete to leading order 
in the coupling constant [§ . For example, the production rate of energetic photons has been 
derived in this way, where a resummed quark propagator containing the effective quark mass 
= g‘^T‘^/6 has been used M- Here we want to reconsider the result of Altherr and 
Ruuskanen ^ performing a complete calculation of the dilepton rate using the Braaten- 
Pisarski method analogously to the photon case. In particular we will study in detail the 
role of m* and M as infrared cutoffs. Like Altherr and Ruuskanen, we will restrict ourselves 
to small photon masses M^T, for which the Os-corrections are of importance, and to large 
photon energies E ^ T. 
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II. DILEPTON PRODUCTION RATE TO ORDER cx'^cxg 


The lowest order matrix elements for the production of virtual photons from the QGP 
are shown in Fig.l. Besides the real contributions (Compton scattering and annihilation) of 
Fig.la,b, appearing also for the production of real photons, there are virtual contributions, 
namely gluon absorption (Fig.lc) and the interference terms of the Born term with radiative 
corrections (self energy insertion and vertex correction) in Fig.Id. Due to phase space 
restrictions and kinematics [E S> T) the gluon absorption as well as the vertex correction 
can be neglected |^. 

The matrix elements can be related to the imaginary part of the photon self energy 
diagrams in Fig. 2 0. Here the processes of Fig.la,b,c can be obtained by cutting through 
the internal quark and gluon lines, while the interference term (Fig. Id) corresponds to a 
cut through the quark lines only. In the latter case the quark self energy becomes on- 
shell, leading to an infrared singularity which cancels the one from the exchanged massless 
quark in the real contribution [Q . This cancellation can be considered as an example of the 
KLN-theorem at hnite temperature |^|. 

The production rate for massless electron and muon pairs derived from the photon self 
energy is given by ^ 

dR 




( 1 ) 


d'^xd'^p bvr^ — 1 “ ^ 

As in the case of the photon production rate M, we want to include the effect of the 
medium by using a resummed quark propagator for soft momenta of the exchanged quark. 


For this purpose we introduce a separation scale kc [0, restricted by gT -C fee -C T in the 
weak coupling limit. For quark momenta smaller than kc we will start from the photon self 
energy shown in Fig.3, where the blob denotes the effective quark propagator, in which the 
quark self energy in the hard thermal loop approximation has been resummed |^. Owing 
to the high energy of the photon and energy-momentum conservation we need to take into 
account only one effective quark propagator and no effective vertices as opposed to the case 
of soft dileptons |T4 . 
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For momenta of the exchanged quark larger than kc, the diagrams of Fig.l or Fig.2, 
containing only bare propagators, will be considered. After adding up the soft and hard 
contributions, the separation scale has to drop out, demonstrating the completeness of the 
calculation |]^. In the case of the photon production a covariant separation scale, i.e. 
lA^I = \uj^ — k‘^\ = kl 1 ^, and a non-covariant one {k = kc) [§[ have been used, both leading 
to the same result. 

In the limit -C ET the hard real contribution (Fig.la,b) of the dilepton production 
rate can be taken over from the photon rate @] , yielding in the case of a non-covariant cutoff 


f dR \ 

y d^xd'^p J 


real 


hard 


10 2 T2 

-r a^as 


, ET 3 In 2 C'(2) 


( 2 ) 


where 7 = 0.57722 and CX2)/C(2) = —0.56996. In the case of a covariant cutoff the constant 
in the square brackets is changed by adding a term In4 — 2 The hard contribution to 
the photon rate follows from a momentum integration over the square of the matrix elements 
of Fig.la,b using distribution functions for the external partons. Applying the Boltzmann 
approximation for the distributions of the incoming partons and using k^. ^ T the result 
(0) is found after some tedious manipulations. It can be obtained more easily by calculating 
the photon absorption rate and applying the principle of detailed balance [|T^ . 

The soft contribution to the dilepton rate can be computed by introducing spectral 
functions for the resummed quark propagator 0 as in the case of the photon rate [ 0 . 
Following the arguments in the calculation of the photon rate we arrive at 


f dR \ 

y d'^xd'^p J 


real 


soft 


40 


d^k 


97r2 M2 J (27r)3 / j dJ 5{E - uj - J) np^uj) 

k) [p+{uj,k)5{uj' + q) + p_{uj,k) 5{uj' - q)] 

+ (1 - q ■ k) [p+(a;, k) - q) + p-(a;, k) + g)]|. 


(3) 


where np is the Fermi-Dirac distribution, q = p — k the momentum of the bare quark 
propagator, and 


4 














p±{uj,k) = 


UJ 




2m*^ 


[(5(a; — u!±) + 6{uj + a;±)] + P±{uj, k) 6{k‘^ — a;^) 


(4) 


the spectral functions with 
^ * 2 


P±{uj,k) = 


ml 


±a; — k 


k(-u,±k)+mf (±l-±g=*lnf±^) 


1 2 


+ 


TT *2 ±u)-k 

2 '"'q k 


1 2 • 


(5) 


The first term of the spectral functions (^), containing the dispersion relations uj±{k) (see 


Fig.4) of the collective quark modes |T^ above the light cone (a; > fc), gives rise to the 
virtual contribution of the dilepton rate. The real soft contribution, on the other hand, 
follows from the second term below the light cone {—k < ut < k), which comes from the 
imaginary part of the hard thermal loop quark self energy. Using E ^ T and assuming uj 
and k to be soft the soft real contribution reduces to 

real ^ 2 


f dR \ 


soft 


^ ^ e y dkdu [{k-uj + E-p) f3+{uj, k) 

+ {k + oj - E + p) , k)]. 


( 6 ) 


The integration range, determined by a; = —k + E — p, oo = k, and the separation scale 
k = kc or k"^ — = k^, respectively, are shown in Fig.4. Compared to the photon case it is 

shifted by U — p > 0, which renders the integration in more difficult. The logarithmic 
dependence on the scale k^ can be extracted by adopting the static limit of the effective 
quark propagator P, 


I3±{uj = 0, fc) 


mf' k 

2 + (7rm*^/2)2 


(7) 


In the case of a non-covariant separation scale kc, where the cn-integration is restricted by 
the limits —k + E — p and A;, while the ^-integration ranges from {E — p)/2 to kc, we find 
to logarithmic approximation 


( dR \ 
\ d^xd'^p ) 


real 


soft 


a 


-EfT 


97r4 M2 


m„ 


In 




+ 4l 


ml 


( 8 ) 


where we have used kc S> gT. In order to determine the function A{E — p) beyond the 
logarithm, we cannot use the static approximation anymore, but have to solve together 
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with d^) numerically. We restrict ourselves to E — p < m*, which holds for realistic values 
of the coupling constant 5 "^!. Then we hnd that A increases from A{E — p = 0) = —1.31, 
which corresponds to the photon result P, to A{E — p = m*) = —1.71. 


Adding the hard contribution (|^) to (||) with = 27ra!sT^/3, the arbitrary separation 


scale kc drops out and we obtain an infrared hnite result for the real contribution of the 
dilepton rate. On the other hand, using a covariant separation scale, (H) is multiplied by a 
factor kl/[{E — pY + k^]. Note that E — p = M ^/{E + p) <^T might be of the same order 
as fee "C T for M T. Hence the separation scale does not cancel in this case in contrast 
to the real photon rate [E = p). This shows that we should demand that uj and k are soft 
individually as it was also assumed in the derivation of the Braaten-Pisarski method, which 
is based on the imaginary time formalism in euclidean space-time 0. 

Finally, we consider the virtual part coming from the pole contribution to the imaginary 
part of the photon self energy in (|l]). Now we have to use a covariant separation scale since 
the infrared singularity in Fig.2 comes from the on-shell self energy at = 0. The soft 
part follows from the pole contribution of the diagrams in Fig.3 corresponding to the hrst 
term of the spectral functions (^. Since uj±{k) lies always below = k^ ^ m*Y (see 
Fig.4), there is no hard contribution to the virtual part. Thus the virtual contribution is 
determined solely by the imaginary part of the self energy in Fig.3 coming from the pole of 
the effective quark propagator, which gives a hnite result when integrated over the entire 
momentum range. 

Combining (§) and the hrst term of we hnd 


( 


d'^xd'^p 


dR 


97r^ M‘^mY P 


5 1 



(9) 


where the limits k^axmin determined from the intersections of the dispersion relations 


'max^min 
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with 10 = ±A; ± {E — p) as shown in Fig.4. Note that the virtnal part va ni shes for M —0 
as k^in tends to inhnity. 

In order to isolate the term proportional to and to investigate its dependence on 
M, m*, E, and T, we introdnce again a separation scale gT -C /c^ -C T. For k < ks we may 
approximate the distribntion fnnctions by ni?(a;±) ~ 1/2 and np{E — ui±) ~ exp(—i?/T), 
whereas for k > kg we may set a;+ ~ fc + m*^/ k leading to o;^ — ~ 2m* ^ and oo- = k [jl^ . 

Hence the plasmino branch uj-{k) does not contribnte for hard momenta as uo- approaches 
k exponentially for k S> gT 0. 

For the part of the integrals in (P) containing the terms with M^/2 it is snflicient to 
restrict to k > kg since the hard part is hnite for —>■ 0. Considering E ^ T the simplih- 
cations k+^^ ~ {E + p)/2, kT^ ~ {E-p)/2, nF{LO+) ~ npik), and nF{E-LO+) ~ npiE-k) 
can be assnmed. Then this term rednces to the Born term, which reads for E ^ T and 
E-p<^T 

/ dR 5 


\ d'^xd^p ) 






( 10 ) 

becanse the effective qnark 


The Born term is contained in the virtnal contribntion ( 
propagator inclndes the bare one. 

Next we consider the terms nnder the integral in (y) proportional to E{uj± ± k). Using 
the approximations discnssed above, we hnd for k < kg 


( dR \ 
y d'^xd'^p ) 


virt 


soft 


a 


-EfT 




dk — k"^) (a;+ — k)+ dk (o;^ — k"^) {tO- + k) 


( 11 ) 


We proceed analogously to Baier et ah in the case of the photon rate, using (o;^ —/c^)(a;±=F 


k)/m* = uj± — k{duj±/dk). Then we obtain with kg S> m* 


f dR \ 
t^d^^xd^p j 


virt 


soft 


Oir-l M2 


14 e-®-"® >71-2 In 


kg 


mo + kz 


+ B 


( 12 ) 


The function B{E — p) decreases from B{E — p = 0) = 0, where the virtual contribution 
vanishes, to B{E — p = m*) = —0.66, where k^^^ = 0 (see Fig.4). 
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The hard part, k > kg, using the approximations given above, reads 


f dR \ 

y d^xd'^p J 


virt 


hard 


10 


*2 

-- 777 

Ovr^ M 2 " 


dk 


nF{k)nF{E — k) 

k ' 


Considering kg E, we hnd 


' dR \ 

^ d^xd'^p J 


virt 


hard 


9M M2 ^ kg 


Adding up both the contributions ([T^) and (p!^ kg cancels and we arrive at 


f dR \ 

i d^xd'^p J 


virt 


a 


-E/T *2 

e ' 


In 


E 


+ kE 


+ B 


( 13 ) 


(14) 


(15) 


contribution vanishes for m* 


Ott^ M2 

We observe that the virtual contribution coming from an interference term is negative, 
thus reducing the dilepton rate compared to the photon rate. Furthermore, the virtual 

0 as it is also the case for the virtual contribution in naive 

H 

perturbation theory [|]. 

The part under the integral in (j^) proportional to (a;± —fc)^ can be shown to be suppressed 
by m*/i? relative to (|I5|) . 

Combining (|^), (|^), and (^) we end up with our hnal result for the dilepton production 
rate 

dR 


d^xd^p 277r^ 


10 2 -E/T i, + k+J ^ 

0 . 0 .g^^ e ' In- - —^-h C 


’M 2 


m: 


(16) 


where the function C ranges from C{E — p = Qt) = —0.73 io C{E — p = m*) = —0.47. 

For giving an approximate analytic expression for kE^ we replace the exact dispersion 
relation ci;+(fc) by a;^ = which deviates from the exact one by less than 11% over 

the entire momentum range. Then we obtain = \Em*^/— M2/(4i7)|, restricted by 


kLn < k+ax - iE + p)/2. 


III. DISCUSSION 

Now we will discuss our result (p!6|) for various limits of M. For M —> 0 the virtual 
contribution vanishes since tends to inhnity for E — p ^ 0. Hence the dilepton rate is 
given by the real contribution which agrees with the photon rate in this limit. 














For M ~ m* we get E — p = M^/ {2E) -C m*. Hence is given by the intersection of 
O' ~ /c + m*^ jk and u) = k + E — p, i.e. {E — p) ~ E. Thus the result for the 

photon rate holds approximately, which agrees with the result of Altherr and Ruuskanen 
1^, if we replace the invariant photon mass in the logarithmic term of their formula (2.18) 
by the effective quark mass in the case of M < m* as suggested in their paper. 


For M r-u T, i.e. E — p ^ T‘^/E, we have to distinguish two cases. First E — p 


m' 


which leads to ~ m* (see Fig.4). Then the logarithm in (|T6|) is given by ln(T/m*). 
Secondly, in the case E — p we hnd from Fig.4 S> m*. Now if ~ we 

recover the photon result, while for ~ T the logarithmic term is given by ln(T^/m*^). 

Finally we will comment on the extrapolation of our result to realistic values of the 
coupling constant, ag = 0.2-0.5, i.e. g = 1.5-2.5. In any case the lowest order result in 
as fails if the rate (|T^) becomes negative. (This unphysical behaviour can originate from 
extracting the leading order contribution and extrapolating to large values of the coupling 
constant. A similar problem has been encountered in the case of transport rates determining 
thermalization times and the viscosity of the QGP BOO Thus one might expect that for 
M ~ T the result gets unphysical, since m* = T for g = a/6. However, for E ':$> T we hnd 
E — p <^T ~ m* and S> T. Hence we end up with a logarithmic term which is always 
positive. However, the constant C behind the logarithm is of the same order (typically 
half of the value of the logarithm), indicating the necessity to go beyond the logarithmic 
approximation. (This is more important for dileptons than for photons due to the negative 
virtual contribution.). 

As an example we have chosen = 0.3 and T = 300 MeV. In Fig.5 the dilepton 
production rate is shown as a function of E for M = 300 MeV and in Fig.6 as a function 
of M for E = 3 GeV. Also shown is the result of Altherr and Ruuskanen [Q. (Here we 
replaced M in their formula (2.18) by m* for M < m*.) For the chosen set of parameters 
both approaches give similar results, since for M < m* = 238 MeV both rates are given by 
the photon result approximately (see Fig.6). However, the dilepton rate found by Altherr 
and Ruuskanen breaks down at somewhat larger values of E as shown in Fig.5. In Fig.5 
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also the Born term O is depicted for comparision. 

Summarizing, we conclude that the complete calculation using a resummed quark prop¬ 
agator according to the Braaten-Pisarski method leads to a hnite result for the dilepton rate 
to order where the effective quark mass m* cuts off the infrared divergence. There is 

no need for using the invariant photon mass M as an infrared cutoff even for M of the order 
of the temperature. 
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FIG. 1. Matrix elements for dilepton production to order a^ag. 
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FIG. 2. Photon self energy diagrams related to the matrix elements of Fig.l by cutting rules. 




FIG. 3. Photon self energy diagrams with an effective quark propagator. 
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FIG. 6. Contribution of order to the dilepton production rate versus the invariant photon 
mass M for T = 300 MeV, a* = 0.3, and E = 3 GeV. 
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